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Results for elastic electron scattering by nuclei, calculated with charge densities of Skyrme forces 
and covariant effective Lagrangians that accurately describe nuclear ground states, are compared 
against experiment in stable isotopes. Dirac partial-wave calculations are performed with an adapted 
version of the elsepa package. Motivated by the fact that studies of electron scattering off exotic 
nuclei are intended in future facilities in the commissioned GSf and RIKEN upgrades, we survey the 
theoretical predictions from neutron-deficient to neutron-rich isotopes in the tin and calcium isotopic 
chains. The charge densities of a covariant interaction that describes the low-energy electromagnetic 
structure of the nucleon within the Lagrangian of the theory are used to this end. The study is 
restricted to medium and heavy mass nuclei because the charge densities are computed in mean 
field approach. Since the experimental analysis of scattering data commonly involves parameterized 
charge densities, as a surrogate exercise for the yet unexplored exotic nuclei, we fit our calculated 
mean field densities with Helm model distributions. This procedure turns out to be helpful to study 
the neutron-number variation of the scattering observables and allows us to identify correlations of 
potential interest among some of these observables within the isotopic chains. 
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I. INTRODUCTION 

Elastic electron-nucleus scattering has been for many 
years a very useful tool to investigate the size and shape 
of stable nuclei [E [3i Hi Electrons interact with 

nuclei basically through the electromagnetic force. If the 
energy of the electrons is high enough, they become a rel- 
atively clean probe to explore precisely the internal struc- 
ture of nuclei, insensitive to strong interaction effects. In 
particular, the analysis of electron scattering data pro- 
vides most valuable information about the charge distri- 
bution in atomic nuclei [1, 0, Hi ■ 

Developments in accelerator technology and detection 
techniques nowadays allow experimentation with nuclei 
beyond the limits of /3-stability. The number of nuclei 
whose masses have been measured keeps growing [9*] and 
this tendency is expected to continue with the use of 
radioactive isotope beams (RIB) [HI [13 • A new 
generation of electron-RIB colliders using stora ge r ings 
is now under construction in RIKEN (Japan) [ij, [lj| 
and GSI (Germany) [l^, [l^. These facilities will of- 
fer unprecedented opportunities to study the structure 
of exotic unstable nuclei through electron scattering in 
the ELISe experiment at FAIR in Germany and the 
SCRIT project in Japan [H Therefore, the theo- 
retical investigation of exotic nuclei with models of pur- 
ported reliability in stable isotopes is a timely and chal- 
lenging problem. Such effort will test the ability of the 
established nuclear theory in the domain of exotic nuclei, 
and may as well provide valuable references for future 
experiments. 

In recent literature, several theoretical studies of elas- 
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tic electron-nucleus scattering in exotic nuclei have been 
reported [13, [U, [H, [H, [H, [2!, [ll] . Some of these works 
are concerned with analyzing electron scattering in light 
nuclei, where exciting exotic phenomena such as the ap- 
pearence of halos may take place. That is the case of 
e^ 6,8jje [11113, "Li [10,[1|, ®B [11[13, and 
[24| nuclei, where possible effects on scattering from the 
occurrence of halos have been investigated. Light nuclei 
require a microscopic treatment of the scattering inter- 
action to pro perly deal with the underlying shell model 
structure [IJ. Other works study the variation of the 
charge form factors along isotopic [13, [ll|, [13| and iso- 
tonic [2^ chains of medium and heavy mass nuclei. It has 
been found that, when the number of neutrons (protons) 
in these isotopic (isotonic) chains increases, the squared 
modulus of the charge form factor and the position of its 
minima show, respectively, an upward trend and a sig- 
nificant inward shifting in the momentum transfer. In 
addition to electron scattering, it is worth mentioning 
that proton scattering may be another valuable tool to 
investigate the changes in the charge density of the nu- 
cleus, especially at its interior, as one proceeds to the 
drip lines along isotopic chains (27j . 

To investigate the internal structure of nuclear charge 
densities, the de Broglie wavelength of the probe has to 
be of the order of 1 fm. This means that the energy 
of the electron beam has to be of the order of several 
hundred MeV. Therefore, for accurate theoretical cal- 
culations of differential cross sections (DCS) and electric 
charge form factors, one needs to solve the elastic scatter- 
ing of Dirac particles in the scalar potential pertaining to 
the nuclear charge distribution. The simplest approach is 
the plane- wave Born approximation (P WBA) , where one 
assumes that the initial and final states of the electron 
can be described by plane waves. Although the PWBA is 
able to account for important features of scattering, it is 



2 



not enough accurate for quantitative calculations of the 
electric charge form factor. A more elaborate method 
is supplied by the Glauber theory using a relativistic 
eikonal approximation of the Dirac equation ^28*1 . It has 
been successfully applied to a systematic study in elastic 
electron- nucleus scattering [23, HE, HH] • The most sophis- 
ticated calculations of electron-nucleus scattering employ 
the exact phase-shift analysis of the Dirac equation. This 
method corresponds to the so-called distorted wave Born 
approximation (DWBA) [2§] and was employed e.g. in 
Refs. [13, Hlj. In the present work we apply a modified 
version of the recently published code elsepa [s^] to the 
elastic electron-nucleus scattering problem. This code 
was originally devised to perform accurate Dirac partial- 
wave calculations of ELastic Scattering of Electrons and 
Positrons by Atoms, positive ions, and molecules in the 
low-energy domain. 

The main input needed for solving the elastic electron- 
nucleus scattering problem is the charge density of the 
target nucleus. In the present article we use charge densi- 
ties obtained with mean field models for calculating elec- 
tron scattering off medium and heavy mass nuclei. For 
this purpose, we employ effective nuclear interactions of 
current use in nuclear structure physics. It is known that 
the overall trends of the elastic electron-nucleus scatter- 
ing in medium and heavy nuclei are, in general, reason- 
ably reproduced by the theoretical charge densities ob- 
tained in the mean field approximation. However, differ- 
ent effective interactions predict electric charge form fac- 
tors and DCS that differ in fine details and describe with 
different quality the experimental data. In our study we 
consider the non-relativistic Skyrme forces SkM* [s^ and 
SLy4 [l^nd the covariant models NL3 [s^, FSUGold 
[13, G2 38], and DD-ME2 They are representative 
examples of effective interactions that accurately describe 
the ground-state properties of finite nuclei and some of 
their collective excitations. 

The parameters of the alluded microscopic interactions 
have been determined from careful calibration to observ- 
ables such as binding energies, single-particle levels, and 
charge and diffraction radii of a variety of selected nuclei. 
The nucleon density distributions of the Skyrme Hartree- 
Fock and relativistic mean field theories are obtained by 
numerical solution of the quantal mean field variational 
equations. The effects of the neutrons on the proton den- 
sity are taken into account, in a self-consistent manner, 
through the interaction terms of the effective force or 
Lagrangian. Therefore, no parametrized shapes of the 
density profiles are implemented. Nevertheless, most of 
the mean field calculations of finite nuclei assume point- 
nucleon densities. The charge density is obtained from 
the proton point-like distribution folded with the proton 
charge form factor [3l| 

PpM = ^e-"^ (1) 

A value c? = 18.29 fm^^ corresponds to a proton root 
mean square (rms) radius of 0.81 fm. This is the standard 



prescription in the fitting procedure of the parameters of 
most Skyrme forces and relativistic mean field interac- 
tions to the experimental data of finite nuclei, including 
charge radii. The covariant models of Ref. 38] arc an 
exception to this fact, and in particular the interaction 
G2 [H, that we will employ in several of our cal- 
culations. The effective Lagrangian of G2 incorporates 
the low-energy electro mag netic structure of the nucleon 
within the theory [H, |48| |. It is to be emphasized that 
in G2 the charge density is obtained directly from the 
solution of the mean field equations and that there is 
no folding to be performed with external single-nucleon 
form factors, hereby maximizing the predictive power. 
In the other mean field forces considered in our work 
we will neglect the contribution of the neutron charge 
form factor [33 | to the charge density. This is known to 
be a reasonabl e ap proximation up to moderate momen- 
tum transfers [33|, which is the region analyzed in the 
present study of scattering. Moreover, it ensures consis- 
tency with the method applied originally to fit the pa- 
rameters of these interactions; additional modifications 
into the charge densities could spoil, e.g., their accurate 
predictions for charge radii. 

It is to be mentioned that the mean field treatment 
would reach its limits in the study of exotic light nuclei, 
where the shell model structure and halos can become 
prevailing features [13]. Therefore we do not attempt to 
treat these light systems in the present work. The devel- 
opment of suitable tools and a unified framework to deal 
with charge densities and electron scattering in nuclei 
across the mass table remains an outstanding problem in 
the field, maybe appropriate for new initiatives like the 
UNEDF collaboration to build a universal nuclear energy 
density functional [ioj . 

The present article is organized as follows. Section 
II is devoted to the theoretical formalism. The mean 
field description of finite nuclei in the non-relativistic 
and relativistic frames is shortly discussed. The Dirac 
partial-wave calculation of elastic scattering of electrons 
by nuclei implemented in the code elsepa is summa- 
rized. In Section III, the elastic electron scattering re- 
sults obtained from the charge densities of the above- 
mentioned effective interactions are studied by compar- 
ing with available experimental data for several stable 
nuclei. In Section IV we investigate the theoretical pre- 
dictions for elastic electron-nucleus scattering in the tin 
and calcium isotopic chains from the proton drip line to 
the neutron drip line. These calculations are performed 
with the relativistic mean field interaction G2 {3^ . Our 
conclusions are laid in the final section. 



II. THEORY 

In the current section we review the basic features of 
the Skyrme and relativistic mean field models that we 
will employ to compute the theoretical nuclear charge 
densities. We also summarize the calculation of Dirac 
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distorted waves for elastic electron-nucleus scattering in 
the code elsepa. The reader conversant with effective 
nuclear mean field models and with knowledge of the 
basics of Dirac partial-wave calculations may prefer to 
proceed directly to the discussion of results that starts in 
Section III. 



A. Mean field description of nuclei 

The mean field approach assumes that nucleons move 
independently in a mean field generated by the other nu- 
cleons of the atomic nucleus. Useful tools for mean field 
calculations of nuclei are the non-relativistic Hartree- 
Fock method with phenomenological interactions and the 
relativistic mean field (RMF) Hartree model with effec- 
tive Lagrangians 41]. These phenomenological forces 
and effective Lagrangians usually depend on about ten 
adjustable parameters that are fitted to reproduce rel- 
evant ground-state properties, such as binding energies 
and charge radii of a few nuclei. 

A common trend of phenomenological interactions 
used in the mean field approach is their simple math- 
ematical structure. In the non-relativistic mean field 
models, the Skyrme interactions [13,113 ^^'^ among those 
most widely used. Skyrme forces are zero-range interac- 
tions, which do not require calculations of exchange con- 
tributions. These forces have been employed for describ- 
ing ground-state properties of the atomic nucleus, low- 
energy excited states, fission and fusion barriers, nucleon- 
nucleus and heavy-ion potentials, etc. (see e.g. Ref. 0). 
The parameter set SkM* is the classical paradigm 
of such a model. It is known to yield charge densities 
in overall agreement with densities inferred from exper- 
iment [l^]. SLy4 [s^ is a more modern version of the 
Skyrme force that was calibrated with special care for the 
isospin sector and for predictions of neutron-rich matter 
that occurs e.g. in neutron stars. 

The RMF theory of hadrons has become another useful 
tool for the study of bulk and single- par ticle properties 
of nuclear matter and finite nuclei p5ll46l. |47| . l48j . In the 
relativistic model, nucleons are treated as Dirac parti- 
cles that interact by exchanging virtual mesons. The co- 
variant theory automatically takes into account the spin- 
orbit force, the finite-range, and the density-dependence 
of the nuclear interaction. The no-sea approximation, 
which disregards effects from the Dirac sea of negative 
energy states, is adopted. The open parameters of the 
model are the meson coupling constants and some of the 
meson masses. After fitting them to binding energies, 
charge radii, and other well-known empirical data of a 
few selected nuclei, the covariant theory predicts aver- 
age properties of spherical and deformed nuclei over the 
whole periodic table in very good agreement with exper- 
iment In, 111,113. 

The original Walecka Lagrangian 45] contained a, oj, 
and p mesons without any meson self-interactions. It 
was able to predict the correct saturation point of nu- 



clear matter, albeit with a very large incompressibility 
modulus. The model was refined with the introduction 
of sigma meson self-interactions (5T| . A parameterization 
of this type is e.g. the celebrated NL3 model [s^]- These 
parameter sets properly describe the data about finite 
nuclei, but often display differences, at densities above 
the saturation point, with microscopic Dirac-Brueckner- 
Hartree-Fock calculations of the nuclear matter equation 
of state [H, [13 • A better agreement with the latter cal- 
culations at densities up to 2-3 times the saturation den- 
sity is achieved by incorporating a quartic vector meson 
self-interaction in the effective Lagrangian. Another in- 
teresting addition is a mixed isoscalar-isovector coupling 
[S^ ]. which allows one to modulate the density depen- 
dence of the nuclear symmetry energy. The variation 
of this coupling leaves the binding energy and proton 
rms of a finite nucleus almost unaltered, but it consider- 
ably modifies the rms radius of the neutron distribution. 
A representative instance of this type of model is the 
FSUGold parameter set t37|. This set yields an equa- 
tion of state that is considerably softer than in NL3, for 
both symmetric matter and neutron matter. Apart from 
the binding energies and charge radii of nuclei, FSUG- 
old delivers a satisfactory description of several modes of 
collective excitations having different neutron-to-proton 
ratios. 

The Lagrangian density associated with the G2 pa- 
rameter set is inspired by effective field theory methods. 
It contains all couplings consistent with the underlying 
QCD symmetries up to the order considered in the ex- 
pansion scheme [H, HI]- In contrast to the majority of 
mean field models, G2 describes the low-energy electro- 
magnetic structure of the nucleon within the theory by 
means of vector-meson dominance and derivative cou- 
plings to the photon, cf. Ref. for details. As indi- 
cated above, this means that no additional calculations 
with external nucleon form factors are needed to obtain 
the charge density, and that the electromagnetic effects 
of the protons and neutrons are included within the low- 
energy regime in a unified framework [s^]- The G2 set 
explains finite nuclei and nuclear matter with a com- 
mendable level of accuracy. It predicts a soft equation of 
state both around saturation and at high densities that is 
consistent with recent measurements of kaon production 
and flow of matter in energetic heavy-ion collisions as 
well as with observations of masses and radii of neutron 
stars [ssj . 

Recent formulations of the RMF theory do not intro- 
duce mesonic self-interactions but make the coupling con- 
stants of the mesons density dependent, like in the DD- 
ME2 parameter set [3§]- These models accurately de- 
scribe the properties of finite nuclei and, in addition, the 
associated equation of state of nuclear and neutron mat- 
ter at supra-saturation agrees with the trends of micro- 
scopic Dirac-Brueckner-Hartree-Fock calculations that 
start from the bare nucleon-nucleon interaction. 

Pairing correlations need to be taken into account for 
the calculation of open-shell nuclei. We will describe 
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them in both, non-relativistic and relativistic frames, 
through a modified BCS approach that takes into account 
the continuum by means of quasi-bound levels owing to 
their centrifugal (neutrons) or centrifugal-plus-Coulomb 
barriers (protons) [56]. For the Skyrme models used in 
this work, the pairing correlations are introduced by us- 
ing a zero-range density-dependent force whose parame- 
ters can be found in Ref . [3^] . In the case of the covariant 
NL3, G2, and FSUGold models, we describe the pairing 
correlations by means of a constant matrix element fitted 
to reproduce the experimental binding energies of some 
selected isotopic and isotonic chains [56] . In the DD-ME2 
parameter set a fixed gap is considered, determined from 
experimental odd-even mass differences (sol ]. 



B. Description of electron scattering 

The ELSEPA code [13] was originally designed for the 
calculation of elastic scattering of electrons and positrons 
by atoms, positive ions, and molecules. We have adapted 
it to handle high-energy electron scattering by nuclei. 
ELSEPA computes the DCS using the conventional rela- 
tivistic partial-wave method, which was first formulated 
by Yennie et al. [1^. The projectile electron is assumed 
to feel the electrostatic field of the nuclear charge distri- 
bution. The potential energy of an electron at a distance 
r from the center of the nucleus is given by 



and 



V(r) = -47re i- 

(2) 

where Pch[r) denotes the charge density of the nucleus, 
considered to be spherically symmetrical. At the energies 
of interest for the electron-nucleus problem, the effect of 
screening by the orbiting atomic electrons is limited to 
scattering angles smaller than 1 degree (see, e.g. Ref. 
[sO]). which are well below the angular range covered 
by the electron-nucleus scattering measurements. Con- 
sequently, electron screening is ignored in the present cal- 
culation. Since the nuclear charge density is assumed to 
vanish beyond a certain radius tb (i.e. the radius of the 
box where the nuclear charge distribution is calculated) , 
the potential ^ is purely Coulombian, V{r) — —Ze^/r, 
beyond that radius. Globally it can be regarded as a 
Coulomb potential with the short-range distortion aris- 
ing from the finite size of the nucleus, i.e., as a modified 
Coulomb potential. 

The DCS for elastic scattering of spin unpolarized elec- 
trons is given by 

1/(0)1' + (3) 
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are the direct and spin-flip scattering amplitudes, respec- 
tively. Here k denotes the wave number of the projectile 
electron, 



chk = x/E{E + 2mec2) 



(6) 



and the functions Pe{cos9) and Pi {cos 6) are Legendre 
polynomials and associated Legendre functions, respec- 
tively. The phase shifts represent the behavior of 
the Dirac spherical waves at large r distances (see e.g. 
Ref. ;5^). 

For modified Coulomb potentials, the spherical solu- 
tions of the Dirac equation are suitably expressed in the 
form 
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The functions f2„^m(f) are the spherical spinors, and the 
radial functions Pek {r) and Q ek [r) satisfy the following 
system of coupled differential equations [53] : 
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dr 
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(8) 



The relativistic quantum number k is defined as k = 
{£ — j){2j + 1), where j and £ are the total and orbital 
angular momentum quantum numbers. Note that j and 
£ are both determined by the value of k; j — \k\ — 1/2, 
£ — j + n/ (21k1). In the numerical calculations, the spher- 
ical waves are normalized so that the upper-component 
radial function Pek {r) oscillates asymptotically with unit 
amplitude. 

For modified Coulomb potentials and r — > oo, we have 
(see, e.g., Ref. (ssf) 



/ TT 

Pek{t') — sin \ kr ^ ^'^ ^ riln2kT 



where 



rj — Ze^mc/{h^k) 



(9) 



(10) 



is the Sommerfeld parameter. It is convenient to express 
the phase shifts 5^ as A, + 5^., where A, is the phase 
shift of the point-nucleus Coulomb potential and 5, is the 
"inner" phase shift of the short-range potential induced 
by the nuclear charge distribution. 

As indicated above, the calculations reported here have 
been performed using the computer code elsepa [30| . 
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It solves the radial Dirac equations using a robust in- 
tegration algorithm, described in Refs. [sl, [s^, which 
effectively minimizes the effect of truncation errors. The 
algorithm starts from a table of values of the function 
rV{r) at the points r,; of a radial grid, which is pro- 
vided by the user. This function is replaced by the nat- 
ural cubic spline that interpolates the tabulated values; 
thus, in the interval between consecutive grid points, the 
potential function rV{r) is represented as a cubic poly- 
nomial. The radial wave equations (O are then solved 
by using the exact power-series expansions of the radial 
functions [5^. The integration is started at r = and 
extended outwards up to a point r^i that is beyond the 
starting radius re of the Coulomb tail. For r > r,n, the 
field is purely Coulombian, and the normalized upper- 
component radial Dirac function can be expressed as 

PE.ir) = cos<5« /^"iW + sinS^g^^ir), (11) 

where fj^^ir) and g'^^^ir) are the upper components of 
the regular and irregular Dirac-Coulomb radial functions 
[S^ . respectively. As usual, the phase shift (5^ is de- 
termined by matching this outer analytical form to the 
inner numerical solution at r^, requiring continuity of 
the radial function Pek{'I') and its derivative. The Dirac- 
Coulomb functions are calculated by using the Fortran 
subroutine described in Ref. ^58'|, which delivers values 
of the regular and irregular Dirac-Coulomb functions and 
their derivatives that are accurate to more than 10 deci- 
mal figures. 

The convergence rate of the series (|4]) and ([5]) for the 
calculation of du/drj is known to be slow. The sum- 
mations are optimized by performing them in two steps 
[30| . First, they are evaluated for the pure Coulomb 
field, for which the phase shifts are known analyt- 
ically and the calculation is fast. Second, the point- 
nucleus results are subtracted from the expansions ^ 
and ([5]); the remaining series represent the effect of only 
the short-range component of the potential and converge 
more rapidly than the original series. As a consequence, 
the number of inner phase shifts one needs to compute 
is normally much smaller than the number of required 
Coulomb phase shifts. 

In general, elsepa gives reliable results for electron 
energies up to 1 GeV but numerical difficulties can ap- 
pear at very small scattering angles because the DCS for 
the bare nucleus is very close to the point nucleus DCS 
(which is known as the Mott DCS) and consequently di- 
verges. The code evaluates the finite-nucleus DCS only 
for scattering angles larger than 10 degrees, which is a 
lower bound of the usually measured electron-nucleus 
DCS. Depending on the particulars of each calculation, 
difficulties can also be found in the high-energy regime 
and for large scattering angles where the DCS takes much 
smaller values than the Mott DCS. In this case, the 
nuclear amplitudes almost cancel the Coulomb scatter- 
ing amplitudes, thus magnifying the numerical errors. 
In practice, round-off errors become apparent when the 



nuclear scattering amplitudes are less then 10~^ times 
the Coulomb amplitudes (which are themselves small). 
When there are indications that these errors could be 
important, elsepa discontinues the calculation and the 
DCS is set to zero. 



III. ELASTIC ELECTRON-NUCLEUS 
SCATTERING IN STABLE NUCLEI 

As stated in the Introduction, we compute the charge 
density distributions with selected effective nuclear in- 
teractions, namely, the Skyrme forces SkM* [s^l and 
SLy4 [ll] , and the covariant parameterizations NL3 [11] , 
FSUGold [13, G2 [H, and DD-ME2 [1|. We obtain 
the DCS from these mean field charge densities using the 
ELSEPA code. Specifically, in the present section we will 
compare the theoretical DCS derived from the mean field 
models with available experimental data for elastic elec- 
tron scattering in at 374.5 MeV [H, 40,42,44,48^^^ 
and ''^Ti at 250 MeV^^'^^Ca at 500 MeV [U, 9°Zr at 
209.6 and 302 MeV ii6,ii8,i24gjj 225 MeV [H, 
and Pb at 248.2 and 502 MeV [H. 

In the conventional analyses of the scattering data 
measured in experiment, the charge density is usually 
modeled by means of an analytical function. For in- 
stance, 2- or 3-parameter Fermi distributions are used 
to this end. The charge density may also be con- 
structed from the eigenfunctions of an adjustable single- 
particle potential of Woods-Saxon or harmonic oscilla- 
tor type. Also, nearly model-independent charge den- 
sities are obtained from a Fourier-Bessel expansion (65j 
with unknown coefficients. In all these cases the free 
parameters in the charge distributions are determined 
from the measured electron scattering data through a 
least-squares minimization procedure. In some compar- 
isons with our theoretical predictions we will employ 
experimental charge densities borrowed from the litera- 
ture. These densities have been extracted from fits with 
Fourier-Bessel expansions d, [3l- In the case of the nu- 
cleus ^^^Sn, for which this type of density is not available, 
we will employ the 2-parameter Fermi distribution given 
in Ref. [63|. The results computed with these "exper- 
imental" (fitted) charge densities will be referred to as 
Exp(fit) in the tables and figures henceforth. 

First, we report in Table [J the rms radii of the theoret- 
ical charge distributions of the nuclei ^^O, *°Ca, ^^Ca, 
^"Zr, ^^^Sn, and ■^'^^Pb predicted by the Skyrme and 
RMF forces used in this work. They are compared with 
the rms radii of the experimental charge distributions ex- 
tracted from the analysis of the electron scattering data 
d, [3, [fill- The agreement is seen to be very good. This 
fact is not surprising because some of the experimental 
values of the charge radius have been used in the fit of the 
free parameters of the Skyrme forces and RMF models 
considered. 

Figure [T] displays the theoretical charge density profiles 
obtained with the investigated mean field models as well 
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TABLE I: Root mean square radii (in fm) for the studied charge distributions. Exp(fit) values are calculated from the 
experimental Fourier-Bessel charge densities @,0|i except for the case of ^^*Sn where a Fermi density is used [63 |. 



Nucleus 


Exp(fit) 


DD-ME2 


G2 


NL3 


FSUGold 


SLy4 


SkM* 


16 Q 


2.74 


2.73 


2.73 


2.73 


2.69 


2.80 


2.81 


«Ca 


3.45 


3.46 


3.46 


3.47 


3.44 


3.51 


3.52 


4«Ca 


3.45 


3.48 


3.45 


3.47 


3.47 


3.54 


3.54 


«OZr 


4.26 


4.28 


4.25 


4.28 


4.26 


4.30 


4.29 


"8Sn 


4.67 


4.63 


4.62 


4.63 


4.63 


4.65 


4.63 


208 


5.50 


5.52 


5.50 


5.52 


5.52 


5.52 


5.51 



G2 

NL3 

FSUGold 

DD-ME2 

SLy4 

SkM* 

. . . Exp(fit) 




r (fm) 

FIG. 1: (Color online) Radial dependence of the charge densi- 
ties of the stable nuclei ^"^O, ^"Zr, and ^"^^Pb. The predictions 
of the different mean field models indicated in the legend are 
compared with the charge densities fitted experimentally by 
the Fourier-Bessel analysis @, 0| • 



as the experimental charge distributions for ^^O, ^''Zr, 
and 208pb gg. In general, theoretical and experimen- 
tal charge densities agree nicely in the fall-ofF region and 
differ more in the nuclear interior as a consequence of 
the shell oscillations of the mean field densities. Differ- 
ences in the inner region are more marked in light nuclei 



where a mean field approximation may be less justified. 
In the present case, if we analyze the region of the ex- 
perimental density of ^^O between the center and about 
2 fm, we see that the covariant sets G2 and NL3 give 
a better description, whereas the SkM* and SLy4 forces 
tend to underestimate it and the DD-ME2 and FSUG- 
old sets tend to overestimate it. A detailed inspection of 
the mean field charge densities shows that they not only 
differ among themselves in the nuclear interior but also 
in the surface region. As discussed in Ref . [3| i the inner 
density is normally larger for nuclei with larger surface 
diffuseness. From Fig. [T] we see that for medium and 
heavy nuclei the Skyrme charge density, in particular the 
one predicted by SLy4, is larger in the interior and con- 
sequently more diffuse at the surface than the RMF dis- 
tributions. Differences regarding the surface diffuseness 
between non-relativistic and relativistic charge densities 
are related to the different density dependence of the ef- 
fective interactions [44| . 

The comparison between theoretical and experimental 
charge densities shall be connected with the discussion 
of the DCS and the electric charge form factors, which 
are the quantities measured in real experiments. Figure 
[5] shows the DCS for elastic electron scattering in ^^O at 
374.5 MeV, ^"Zr at 302 MeV, and ^ospb at 502 MeV, 
which we choose as representative examples to illustrate 
the predictions of the different mean field models. To 
improve readability, the curves calculated with NL3 and 
SkM* are not displayed in the figure. In the three scatter- 
ing processes, one observes that the experimental data, 
are similarly reproduced by all of the considered mean 
field charge densities at small scattering angles, up to the 
first diffraction minimum. The statement is more valid 
for lead, where models and data keep close up to the sec- 
ond, or even third, diffraction minimum. The case which 
is seen to pose more difficulties to the mean field models 
at all scattering angles is, not surprisingly, the lightest 
investigated nucleus, ^^O. From the considered interac- 
tions, only the parameter set G2 is able to reproduce the 
experimental data of ^^O at 374.5 MeV closely at all scat- 
tering angles. As one might expect from the discussion of 
Fig. [U the deviations among the DCS predicted by the 
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FIG. 2: (Color online) Elastic DCS for electron-nucleus scat- 
tering in ^'^O, ^"Zr, and ^°*Pb as a function of the scattering 
angle 8, at the energies shown. The results from the mean field 
models indicated in the legend are compared with the mea- 
sured DCS (Exp) [63,[63,[6a| and with the DCS calculated in 
the ELSEPA code from the charge densities fitted experimen- 
tally by the Fourier-Bessel analysis [Exp(fit)] 0, 0|- 




FIG. 3: (Color online) Squared charge form factor for ^^O, 
^''Zr, and ^'^^Pb as a function of the momentum transfer q, at 
the electron beam energies shown. It has been obtained by 
applying Eq. psp as described in the text, both for the mean 
field models indicated in the legend and for the experimentally 
fitted charge densities [Exp(fit)] [1, 



various mean field models, and the discrepancies with 
respect to experiment, become more prominent at the 
largest scattering angles in the three processes studied in 
Fig. [21 One detects a significant difference between the 
non-relativistic Skyrme interactions and the RMF mod- 
els used here. It is seen that the SLy4 interaction yields 
DCS values that, in general, are smaller than those cal- 
culated using the charge densities of the covariant sets 
(the same happens with the Skyrme force SkM*). This 
trend is especially clear in ^^O and ^°Zr, if one just ex- 
cludes the region immediately after the first diffraction 
minimum where a crossing of the DCS values predicted 
by some of the models takes place. 

One of the most effective quantities to characterize 
elastic electron-nucleus scattering is the electric charge 
form factor F(q). The momentum transfer q is related 
to the scattering angle 9 in the laboratory frame by 

chq = 2E sm{e/2). (12) 

In PWBA the electric-charge form factor is computed 
as the Fourier transform of the charge density. In the 
present work we obtain |i^((7)p at a given beam energy 
as 



where da/dO. is the DCS calculated from the DWBA 
analysis and dayi/dVl is the Mott DCS. Equation (fTS]) 
goes beyond the first Born approximation to the charge 
form factor because, instead of the PWBA for the DCS 
of the point nucleus, the exact Mott DCS is used. It 
will be referred to as -FbwBA('?) in what follows. We will 
extract the squared charge form factor from experiment 
using the same expression (|13l) . by inserting the DCS ob- 
tained from the experimentally fitted charge densities on 
its right-hand side. 

Figure [3] displays the g-dependence of |^bwBA(9)|^ for 
the nuclei ^^O, ^"Zr, and ^°®Pb at the electron beam 
energies 374.5, 302, and 502 MeV, respectively. Results 
are calculated from the considered mean field models and 
from the experimental charge density (as in Fig. [51 the 
curves obtained with NL3 and SkM* are not shown) . As 
it may be expected from our previous analysis of the 
DCS, the experimental values of |F((7)|^ are well repro- 
duced in the low-momentum transfer region by all of the 
discussed mean fields. However, some discrepancies ap- 
pear between the theoretical predictions and the exper- 
imental data at large momentum transfers (or, equiva- 
lently, at large scattering angles). They are an indication 
that the theoretical mean field models describe differe ntly 
the central region of the experimental charge density [60| , 
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FIG. 4: (Color online) Relative differences of DCS in pairs of neighbor nuclei at 250 MeV for Ca and Ti and at 225 MeV for 
Sn. The results from mean field models are compared with the measured RDDCS (Exp) and with the RDDCS calculated in 
the ELSEPA code from the charge densities fitted experimentally [Exp(fit)] in Refs. [6ll. l63l]. Note that in the present figure the 
vertical scales are linear and that all of them have been magnified by a factor 100. 



which we have addressed in the discussions of Fig. [TJ 

A more exigent test of the mean field charge densities, 
than the DCS themselves, is provided by the analysis of 
the relative differences of DCS values between pairs of 
neighbor nuclei (RDDCS). The RDDCS are defined as 



D{A - B) 



{da/dn)A + {da/dn)B' 



(14) 



We explore the following cases where experimental val- 
ues are available: ^^Ca-'^'^Ca, 4°Ca-'**Ca, ''°Ca-'*sCa, and 
48Ca-48Ti [m, as well as "^gn-HSSn and ^^^Siv-^'^^Sn 
[ill. The theoretical predictions for the RDDCS ^ 
from the mean field charge densities are displayed against 
experiment in Fig. 01 where a factor of 100 has been in- 



troduced for the sake of clarity. Notice that the vertical 
scale of this figure is linear instead of logarithmic. In 
general, all the considered models describe fairly well the 
experimental RDDCS values for small scattering angles. 
The agreement with experiment deteriorates when the 
scattering angle increases, pointing out again some pos- 
sible deficiencies in the inner region of the theoretical 
charge density distributions. 

To make a more quantitative analysis of the electron 
scattering DCS derived from the mean field charge den- 
sities, we introduce a normalized square weighted differ- 
ence (rf^), or comparison magnitude, with respect to the 
DCS measured in experiment. For each nucleus and elec- 
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TABLE II; Normalized square weighted difference (d^) between the calculated and measured DCS values. Exp(fit) values are 
calculated from the experimental Fourier-Bessel charge densities @, except for the case of ^^*Sn where a Fermi density is 
used ■ The energy E of the incident electrons is in MeV. 



Nucleus E Exp(fit) DD-ME2 G2 NL3 FSUGold SLy4 SkM* 



16(3 


374.5 


11.1 


«Ca 


250.0 


7.18 




500.0 


3.48 


«Ca 


250.0 


6.66 




500.0 


3.19 


9»Zr 


209.6 


0.78 




302.0 


0.86 


"«Sn 


225.0 


5.43 




248.2 


30.6 




502.0 


21.2 



88.7 


13.1 


38.6 


3.15 


16.2 


13.9 


1.49 


42.9 


19.7 


4.85 


9.74 


7.14 


1.11 


17.0 


3.53 


0.87 


2.21 


1.36 


0.91 


9.92 


3.27 


18.4 


34.8 


25.5 


44.4 


154. 


74.8 


14.1 


186. 


50.5 



206. 


191. 


194. 


0.84 


24.4 


24.3 


5.79 


40.0 


39.0 


4.08 


14.9 


13.6 


2.57 


21.84 


18.5 


0.65 


6.53 


5.36 


0.67 


9.35 


7.19 


31.8 


2.75 


4.20 


89.5 


89.2 


61.0 


61.1 


95.9 


76.5 



TABLE III: Normalized square weighted difference (d^) between calculated and measured relative differences of DCS in pairs 
of neighbor nuclei. The beam energy per electron is 250 MeV for the Ca isotopes and **Ti, and 225 MeV for the Sn isotopes. 



D{A-B) Exp(fit) DD-ME2 G2 NL3 FSUGold SLy4 SkM* 



D(«Ca-*2Ca) 0.56 

D(«Ca-*''Ca) 1.14 

D(«Ca-**Ca) 1.06 

75(48ca-^»Ti) 2.49 

D("^Sn-"**Sn) 2.05 

D("«Sn-"''Sn) 4.03 



9.1 


28.3 


16.0 


4.5 


29.6 


12.2 


16.4 


4.89 


7.74 


18.0 


19.6 


31.0 


8.05 


7.80 


9.00 


5.35 


6.98 


7.50 



11.1 


9.11 


12.9 


3.88 


7.08 


9.13 


38.5 


94.1 


49.3 


37.8 


71.8 


64.9 


10.1 


13.2 


18.5 


9.22 


7.05 


7.18 



tron beam energy, it is defined as 



1 ^ 

1=1 



(dcr/dr2)f P - {da/dn)1 
5{da/d^lf^' 



calc 



(15) 



In this expression, the quantities [da / dU,)1^'^'^ , {d<j/dfl)1^^ 
and 6{da/dD,)'^^^ are, respectively, the calculated DCS, 
the measured DCS and the uncertainty of the latter. The 
sum in (fT5|) runs over all the N available data for the 
given scattering process. 

In Table [Til we report, for several elastic electron- 
nucleus reactions, the values (fTS]) from theory in com- 
parison with the obtained from the charge densities 
fitted experimentally. One realizes that the of the 
DCS computed with the theoretical charge densities take 
sizably varying values, depending on the scattering pro- 
cess and nuclear model. Nevertheless, for all the con- 
sidered electron scattering processes, there are Skyrme 
forces or RMF parameterizations that yield a d^ of sim- 
ilar quality to the experimental charge density. In par- 
ticular, the RMF sets DD-ME2 and FSUGold give, on 
the average, the best overall agreement with the consid- 



ered experimental data. The Skyrme forces SkM* and 
SLy4 provide a similar description, but globally this de- 
scription is slightly worse than in the case of the RMF 
parameterizations, except for the nucleus ^^^Sn. 



We have also analyzed the d^ values of the relative dif- 
ferences between differential cross sections in pairs of 
neighbor nuclei. The corresponding results are collected 
in Table lllll For all the mean fields considered in the 
present study, the predicted d^ values are larger than 
the values computed with the charge densities that 
are fitted to experiment. Looking globally at the results 
presented in Table IIIIl the agreement obtained between 
theoretical and experimental RDDCS is not significantly 
worse than the agreement found in the case of the DCS 
analyzed previously in Tabic HIl But, in the DCS analy- 
sis, it has been seen in Table [ll] that for each scattering 
process there is always one or various specific interactions 
that are able to come very close to experiment, whereas 
in Table IIIIl this does not happen in any of the studied 
RDDCS in pairs of neighbor nuclei. 
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IV. ELECTRON SCATTERING ALONG 
ISOTOPIC CHAINS 

Different tendencies of the square of the charge form 
factor as a function of the momentum transfer have been 
studied in earher hterature for isotopic [13, HH, [2^ and 
isotonic [26l | chains of medium and heavy nuclei. Our 
aim here is to perform an analysis, quantitative when- 
ever possible, which could eventually be useful for future 
electron scattering measurements in RIB facilities. 

We will be concerned with the study of electron scat- 
tering in isotopic chains of medium and relatively heavy 
systems, which will be exemplified by the cases of the Ca 
and Sn isotopes. Many of these nuclei lie in the region 
of the nuclear chart that is likely to be explored employ- 
ing RIB facilities. Some of them may be investigated 
in_future electron scattering experiments, such as ELISe 
^ and SCRIT [H, [l|. For the purpose of our study, 
we are interested in predicting the trends of the varia- 
tion along the isotopic chains of the electric charge form 
factor in the low-momentum transfer regime. Our calcu- 
lated mean field charge densities will be parameterized by 
means of the Helm model [g^ , often used in the analysis 
of experimental data, with a view to gain deeper physi- 
cal insights and to elucidate possible correlations among 
scattering observables within the isotopic chains. 

We have seen in the previous section that in the re- 
gion of small q values the experimental results are almost 
equally well reproduced by the calculations with the dif- 
ferent theoretical mean field models considered. In our 
subsequent study of electron scattering in the tin and cal- 
cium chains, as a representative reference, we will work 
with the charge densities predicted by the covariant in- 
teraction G2 [sl, . This model was constructed as an 
effective hadronic Lagrangian consistent with the sym- 
metries of quantum chromodynamics. As mentioned, the 
model describes the low-energy electromagnetic structure 
of the nucleon using vector-meson dominance and pro- 
vides directly the charge density of the nucleus, so that 
no external single-nucleon form factors are required to 
compute the latter G2 is also a rehable parameter 
set for calculations of ground-states of nuclei and, at the 
same time, for predictions of the nuclear equation of state 
up to supra-normal densities and of some properties of 
neutron stars [H, HI, [5^1 ■ Calculations of the squared 
charge form factor done in PWBA with the set G2 for 
stable isotopes have been reported elsewhere [38]. 

For each nucleus in an isotopic chain we compute the 
associated DCS via the DWBA calculation using the G2 
charge density. The electron beam energy in all the in- 
vestigated scattering processes is fixed at 500 MeV. In 
practice, the energy dependence of the electric charge 
form factor defined in Eq. is seen to be consider- 
ably weak for low-momentum transfers, the regime ad- 
dressed in our analysis. We have verified numerically 
that |-FbwBA(9)P of high-energy electron scattering de- 
pends little on the electron beam energy for momentum 
transfers up to 1-1.5 fm~^ (the precise value depends on 




q (fm ') 

FIG. 5: (Color online) Squared charge form factor of ^^*Sn 
and ^°*Pb, derived from the covariant mean field model G2, 
for the indicated electron beam energies. The figure points 
out that the sensitivity of the charge form factor extracted 
through Eq. (|13|) to the beam energy in high-energy elas- 
tic electron-nucleus scattering is rather small, even when the 
atomic number of the target is large as in ^"^Pb. 



the nucleus). This happens even for a heavy system like 
^•^^Pb, where the departure from the point-nucleus as- 
sumption of the Mott DCS used on the r.h.s. of Eq. ([T3|) 
is more significant. We illustrate the situation in Fig. [5l 
where we display |-FbwBA(9)|^ at several beam energies 
for electron scattering off ^^^Sn and ^°*Pb. 

Before proceeding to the presentation of the results in 
the tin and calcium isotopic chains, we briefly summarize 
in the next subsection how we determine the parameters 
of the Helm model charge density distributions. 

A. Equivalent Helm charge densities 

Valuable insights into the study of electron scattering 
often stem from consideration of modeled charge densi- 
ties and electric charge form factors. Moreover, these pa- 
rameterized forms are instrumental in experimental data 
analyses. A notable case is the so-called Helm model, 
whose original version ^6d\ has later been extended in 
various ways for more accurate descriptions of the exper- 
imental charge densities [g^, [H, [6^ . In the simpler ver- 
sion of the model, two chief features of the nuclear charge 
density, namely, the position and the thickness of the sur- 
face, can be related explicitly to the electric charge form 
factor obtained in PWBA. The Helm charge density is 
obtained from the convolution of a constant density po in 
a hard sphere of radius Rq (the diffraction radius) with a 
Gaussian distribution of variance (whose square root 
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relates to the nuclear surface thickness): 

pW(f) = / dr'/G(f - r')poe(i?o - r), (16) 



where 



(17) 



The rms radius of the Helm density is readily obtained 
from Eqs. (fT6|) and (fT7|) . It can be expressed in terms of 
Rq and a as 



(18) 



The corresponding electric charge form factor in PWBA 
is given by 



qRo 



(19) 

where ji(x) is a spherical Bessel function. 

The diffraction radius Rq of the Helm density is usu- 
ally fixed as follows. One requires that the first zero of 
Eq. (fT9|) occurs at QminRo , where qmin corresponds to the 
first minimum of the modulus of the PWBA form factor 
(FpwBA (q) hereinafter) associated to the original charge 
distribution that the Helm density attempts to describe: 



Rq 



4.49341 

^min 



(20) 



The variance cr^ of the Gaussian is chosen to reproduce 
the height of the second maximum of |fpwBA('?)|7 located 

"it Qmax- 



■In 



3jl(9max^o) 



n 



(21) 



For moderate values of qRo, the Helm charge form fac- 
tor i^(^) (g) reproduces well the actual charge form factor, 
with the exception of the regions closest to its zeros . 
The relative difference between F^^\q) and the actual 
charge form factor becomes progressively manifest as the 
momentum transfer grows. The applicability of the Helm 
model near the drip lines is to be explored. The extent 
to which it may be appropriate away from stability, for 
the purposes of our study, will be validated later from the 
numerical point of view in connection with the discussion 
of the calculations in the Sn and Ca chains. 



B. Tin and calcium isotopic chains 

We are interested in the study of elastic electron scat- 
tering in the Ca and Sn chains. The calculated DCS and 
squared charge form factors |i^DWBA(9)P show, for the 
lightest nuclei considered here (calcium isotopes), a rel- 
atively well marked first minimum. This first minimum. 



however, practically disappears for the heavier nuclei an- 
alyzed (tin isotopes). (This fact can also be told from 
the previous figures [H and E] for i^Q, 90Zr, and ^ospb.) 
In the latter case, the form factor |-FbwBA('?)P of Eq. (fT5|) 
still shows an inflection point (IP) in the low- momentum 
transfer region, a point where the curvature changes sign. 
In the absence of an explicit minimum at low-momentum 
transfer, these IP are the best candidates to characterize 
along the isotopic chain relevant properties of the electric 
charge form factor in the small-g region. 

We determine for each isotope an equivalent Helm dis- 
tribution from the calculated mean field charge density. 
In the PWBA the distortion of the electron wave func- 
tions due to the Coulomb potential of the nucleus is ne- 
glected. The effect of the Coulomb attraction felt by the 
electrons can be simulated by replacing the momentum 
transfer q by an effective value (66[ 



9cff = g 1 



Za 



qR, 



ch 



(22) 



where Rch — \/oj5R is the rms of the charge density as- 
suming a hard sphere distribution of radius R = tqA^^^. 
In Ref. I&S] the value of the constant c is taken to be 3/2. 
Here we leave c as a free parameter. It is optimized so 
that the rms radii of the equivalent Helm charge densities 
P^ . with Rq and a determined by Eqs. ^ and (PT|l . 
and with q replaced by qcs , best reproduce along the iso- 
topic chain the mean field rms charge radii obtained with 
the RMF parameterization G2. Proceeding in this way, 
we find c « 0.15 for the tin isotopes and c « 0.12 for the 
calcium isotopes. 

We first investigate the tin isotopic chain. The calcu- 
lated Helm parameters i?o and ct^ are displayed in the 
lower and upper panels of Figure [SI respectively. It is 
seen that Rq steadily increases with the mass number 
A and that it roughly follows the typical A^^^ law. On 
the contrary, the change of the variance cr^ with mass 
number shows a non-uniform character along the isotopic 
chain, related to the underlying shell structure. It oscil- 
lates around a mean value ~ 0.72 {a ~ 0.85). It is 
to be noted that displays local minima at the dou- 
bly magic isotopes ^'^^Sn and ^^^Sn (neutron drip line 
nucleus). This fact points out a stiffer nuclear surface 
for neutron-rich nuclei with double-closed major shells, 
in agreement with earlier literature [67| . 

The upper panel of Figure [7] depicts the radial depen- 
dence of the mean field charge density profiles computed 
with the G2 covariant interaction for the isotopes ^°"Sn, 
^^^Sn, and ^^^Sn. We have selected these examples to 
illustrate the evolution of the results along the isotopic 
chain from one drip line to another. Though the three 
isotopes share the same atomic number, one notices out- 
standing variations in the calculated mean field charge 
densities of these isotopes, both in the surface region 
and, especially, in the interior region. This fact reflects 
the important influence of the changing neutron number 
on driving the structure of the charge density along the 



12 



0.75 

I, "•■^ 
0.65 



<o"> = 0.72 fm 



100 110 120 130 140 150 160 170 
A 



6.4 
6.2 

I ^ 
5.8 

5.6 

5.4 



Tin Isotopes 



0.08 



0.06 = -J^ 



tl. 0.04 



0.02 - 




5.2 



4.8 



4.6 



4.6 



4.8 



5.2 



5.4 



5.6 



o 


'1 

G2 


' 1 ' 1 


X 


Helm 


















M K 


Tin Isotopes 




,1 


,1,1 


6 


4.8 5 


5.2 5.4 5.6 



FIG. 6: Mass-number dependence of the Helm parameters 
(upper panel) and i?o (lower panel) predicted by the covariant 
mean field model G2 in the Sn isotopic chain. The average 
value of is depicted by a horizontal dashed line. 



FIG. 7: Charge densities (upper panel) and charge radii 
(lower panel) in the Sn isotopic chain, according to the co- 
variant model G2 and to the determined Helm distributions. 



isotopic chain. This influence is encoded in the inter- 
action terms of the covariant Lagrangian of G2 through 
the exchanged mesons and the couplings to the photon, 
and it is accounted for self-consistently by the mean field 
calculations. In the same figure we compare the fitted 
Helm charge densities and the original G2 charge den- 
sities, for the three discussed tin isotopes. At the in- 
terior of these nuclei, we see that the uniform density 
of the Helm model averages the oscillations of the self- 
consistent quantal densities obtained with the G2 inter- 
action. In spite of the fact that the surface fall-off of the 
Helm densities is of Gaussian type, the agreement at the 
surface region between the mean field charge densities 
and their Helm equivalents is fairly good as one proceeds 
along the isotopic chain, from stability to the neutron 
and proton drip lines. This provides some confidence on 
using the Helm model as determined in the present work 
when the drip lines are approached in the Sn chain. The 
lower panel of Figure [7] illustrates the mass number de- 
pendence of the mean field and the equivalent Helm rms 
charge radii. Our approach leads to an excellent agree- 
ment between the Helm values and the mean field values 
in the Sn chain. One appreciates some slight discrepan- 
cies only for isotopes very close to the drip lines. The rms 
charge radii of both calculations follow the expected lin- 
ear trend with A^^^ . One observes some departure from 
the A^/'^ behavior to slightly lower values in the isotopes 
close to ^ = 132. This is finally the reason why the vari- 
ance (7^ decreases around A = 132 in Figure [S] [also see 



Eq. (HID]. 

We now analyze the evolution along the isotopic chain 
of the momentum transfer qip corresponding to the first 
inflection point of |i^DWBA(<z)|^ and the variation of 
-FbwBA('?ip)P, calculated for an electron beam energy of 
500 MeV. We discuss possible correlations with the pa- 
rameters i?o and (T of the equivalent Helm density. Two 
noticeable findings of this study are displayed in Fig. [S] 
In the upper panel of the figure it is seen that the change 
of the momentum transfer qip with mass number in the 
Sn isotopic chain approximately follows an A~^/^ law. 
Thus, the position of the first minimum of the square of 
the modulus of the electric charge form factor is shifted 
towards smaller q values as the number of neutrons in the 
isotopic chain increases, a feature also noted in Ref. [25l |. 
Another interesting result is that the momentum transfer 
qiP turns out to be almost proportional to the effective 
momentum transfer q^s corresponding to the corrected 
first minimum of |^pwba('?)| computed with the mean 
field charge density. The situation is illustrated in the 
middle panel of Figure [H Recalling Eq. ((20)) . this corre- 
lation allows one to establish a straightforward relation- 
ship between gip and the Helm parameter i?o in the tin 
isotopic chain: 

4 934 

i?o « (23) 
gip 

We also find that the rms charge radii of the tin isotopes 
exhibit a considerable linear correlation with the momen- 
tum transfer gip , as depicted in the lower panel of Fig. [S] 



13 



0.9 



0.85 



7.4^ 
7.2 - 

7 - 



r = 0.992 



0.8 - 



0.18 



0.19 



0.2 



0.21 



0.84 

'a 0.8 

I 0.76 



0.72 
5.2 

5 

4.8 

4.6 

4.4 



' 1 ' 1 ' 


1 1 1 1 , 


/■ = 0.997 






Tin Isotopes 


] > 1 > h 

» i 


— 1 \ \ 1 ^ 

G2 
X Helm 




« o " 

^ 



0.78 0.8 0.82 0.84 0.86 0.88 0.9 
?,p (fm~') 



^- 6.8 - 
I 6.6- 
3 6.4 - 
6.2 - 
6- 



Tin Isotopes 
500 MeV 



0.4 



0.45 



0.5 

2 2 



0.55 



0.6 



FIG. 9: Correlation along the Sn isotopic chain of the squared 
charge form factor in DWBA, at its first infiection point (gip), 
with the product of the Helm parameter and gjp . The re- 
sults are computed with the G2 effective interaction for scat- 
tering by 500 MeV electrons. Note that each point of the 
figure corresponds to a different isotope and that the mass 
number increases from right to left. 



FIG. 8: Results predicted by the G2 effective interaction for 
500 MeV electron scattering in the Sn isotopic chain. Upper 
panel: mass-number dependence of the momentum transfer at 
the first inflection point (qip) of the squared charge form fac- 
tor in DWBA. Middle panel: correlation of the effective mo- 
mentum transfer at the first minimum of the squared charge 
form factor in PWBA (geff.min) with the value of qip. A linear 
fit of the results is shown and the correlation coefficient r is 
indicated. Lower panel: the change of the charge radii calcu- 
lated with G2 and with the corresponding Ifelm densities is 
depicted against the value of gip. 



These correlations could, in principle, provide an alter- 
native way to obtain the parameter a of the equivalent 
Helm charge densities directly through Eq. p8|) . taking 
into account the relationship (|23p between Rq and qip. 

Figure [5] displays the evolution in the Sn chain of 
|^bwBA(9ip)P (the value of the DWBA squared charge 
form factor calculated at qip) as a function of a^qfp- The 
addition of neutrons along the isotopic chain in general 
brings about an increase of the value of |fDWBA(qip)P, 
also documented in the literature [1^ [2l], [51, |24|, l2a |. 
Furthermore, we notice that an interesting linear corre- 
lation arises between the quantities |i^DWBA(9ip)P and 
CT^QiP as one moves from the proton-rich side of the iso- 
topic chain to the neutron-rich side. This correlation may 
be qualitatively understood in the following terms. If for 
guidance we consider the expression (|19p of the electric 
charge form factor in the adopted Helm model, we see 
that the natural variables to investigate the variation of 
the charge form factor are qRo and a^q^. But, as stated 
in Eq. the value of qRo at the first inflection point 
of |-FbwBA(<z)|^ is practically independent of the mass 



number in the whole isotopic chain. Therefore, the mass- 
number variation of cr^^jp is left as the principal source 
for the ^-dependence of the value of the squared charge 
form factor at q = qip. It is then reasonable that the 
change with A of |Fdwba(9ip)P and (J^qfp is correlated 
along the isotopic chain, a feature confirmed by Fig.[9l As 
a physical insight, helped by Eq. (^5]) . the product (J^qfp 
can be recast as proportional to /R^, which is the ratio 
between the surface width and the mean location of the 
surface of the underlying nuclear charge density. 

Similar results to those discussed above for tin have 
been obtained in the study of the calcium isotopic chain. 
They are presented in Figures [TUl [T51 We first have mod- 
eled the calculated mean field densities of the calcium 
chain by the simpler form of Helm densities. The mass- 
number evolution of the parameters of the Helm charge 
densities is illustrated in Figure [TOl The shrinkage of the 
Helm parameter a at magic neutron numbers, noticed 
in the tin chain, is a feature also present in the calcium 
isotopes for the magic neutron numbers N = 28 and 
50. However, the effect at = 20 is completely washed 
out, in agreement with the result of Ref. [63|- The mean 
field charge densities of ^^Ca, ^^Ca, and ™Ca calculated 
with G2, along with the equivalent Helm profiles, are dis- 
played in the upper panel of Figure 1111 As in the case of 
the tin chain, the effects of the addition of neutrons are 
very manifest in the mean field charge density, as one can 
appreciate from the changes in the interior and surface 
regions of the density distributions of the shown Ca iso- 
topes. One also sees that the Helm densities manage to 
follow on average these changes. The Helm profiles are 
found to reproduce closely the surface region of the mean 
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FIG. 10: Same as Fig. |6] for the Ca isotopic chain. 



FIG. 11: Same as Fig. [7] for the Ca isotopic chain. 



field densities, including the isotopes at the proton- and 
neutron-rich sides of the calcium chain. Discrepancies 
are observed in the interior region of the density distri- 
butions. In the lower panel of Figure [Tl] we display, as 
a function of A^^^, the variation in the calcium chain of 
the rms radii of the mean field and Helm charge den- 
sities. Compared to the case of Sn, this lighter chain 
presents more significant departures from the A^^^ be- 
havior. Also, the agreement of the rms radii of the Helm 
model with the self-consistent G2 values is less good. 
This is more visible in approaching the drip lines, es- 
pecially at the proton drip line. 

Analogous correlations to those previously discussed in 
the case of tin, between the value of q at the first mini- 
mum (or IP) of |^bwBA(9)P with (i) the mass number, 
(ii) the first minimum of |i^pwBA(Q)|, and (iii) the rms 
radius of the charge distributions, are similarly found in 
the analysis of the calcium isotopes. They are displayed 
in Figure [T^ In the present case, however, a significant 
departure of the radii i?o and i?ch from the A^^^ law 
is observed as one moves towards the proton drip line. 
These deviations may be largely due to the fact that, in 
approaching the proton drip line, the protons are more 
loosely bound and therefore the charge density extends 
to larger distances compared with the stable nuclei above 
^"^Ca. The effect is much more prominent in calcium than 
in tin because of the lower Coulomb barrier due to its 
smaller atomic number. In turn, the same effect may 
originate that the Helm parameter a of ^°Ca, measuring 
the surface thickness of the nucleus, does not decrease 
as compared with the heavier neighbor nuclei (see upper 



panel of Figure [TU)) . We plot in Fig. [T3] the square of 
the DWBA charge form factor computed for 500 MeV 
electrons at its first minimum against the value of the 
product cr^g^jjj. In calcium, as in the case of the tin iso- 
topes, an outstandingly linear correlation exists between 
both quantities. 

Finally, in order to validate the consistency of our anal- 
ysis with Helm density equivalents, we compute in test 
cases the squared DWBA electric charge form factor both 
with the Helm profiles and with their original mean field 
charge densities. The results, as a function of the mo- 
mentum transfer g, are compared in Fig.[T3]in a few iso- 
topes of calcium and tin. The cases shown in the figure 
are chosen to illustrate the situation from the proton to 
the neutron drip lines, but similar conclusions are found 
for the other isotopes of these chains. In the range of q 
values up to Ril.5 fm~^ in the calcium isotopes and up 
to around ~1 fm^^ in the tin isotopes, one observes an 
excellent agreement between the results for |Fdwba(<7)P 
obtained using the mean field charge densities and using 
the equivalent Helm charge densities. A similar situa- 
tion is found in the nuclei that lie near the drip lines: the 
agreement between the Helm and the mean field results is 
slightly worse but without sizable differences compared to 
the more stable nuclei. This scenario points out that the 
electron-nucleus scattering results in the low-momentum 
transfer region computed with mean field charge densi- 
ties, are very well simulated by equivalent Helm densities 
with the Rq and a parameters obtained according to the 
method described above. As the value of the momentum 
transfer increases, expectably, the discrepancies between 
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the actual calculations with G2 and the adjusted Helm 
model become evident. 

Let us now concentrate on the results pertaining to 
the self-consistent calculations with G2 that are shown 
in Fig. 1141 One can observe that, as we have mentioned 
earlier in the text, the first minimum that is visible in 
the squared charge form factor of calcium sinks away in 
tin. Its fingerprint in tin is recognized as an inflection 
point. The second and further minima of |^bwBA(9)P 
remain, however, clearly identifiable in the tin isotopes. 




q (fm" ) 

FIG. 14: (Color online) Comparison in Ca and Sn isotopes 
between the squared charge form factor in DWBA calculated 
with the charge densities of the G2 mean field interaction 
(solid lines) and with the equivalent Helm charge densities 
(dashed lines). Notice that the vertical and horizontal scales 
for calcium and tin are different. 



With the progressive addition of neutrons in either of the 
two isotopic chains, it may be seen from Fig. [14] that the 
first minimum of |-FbwBA('?)P, or its signature, becomes 
slightly more marked. This effect is, however, much less 
noticeable than the discussed effect on the first minimum 
induced by changing Z. One also observes that the lo- 
cation of the minima or inflection points of the squared 
charge form factor is gradually pushed towards lower mo- 
mentum transfers as the nucleus becomes more and more 
neutron rich. This effect has been noted in previous lit- 
erature [m, m, [13, m, and for the first minimum 
we have described it in more detail in the upper pan- 
els of Figs. [5| and [T^l The same trend (inwards shift of 
the minima or IP) occurs with increasing atomic num- 
ber (compare the scales of the q axis for Ca and Sn in 
Fig. In general, the inwards displacement of the 

momentum transfer of the minima is accompanied by a 
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simultaneous increase of the height of the maxima of the 
squared charge form factor. 



V. SUMMARY AND CONCLUSIONS 

We have adapted the elsepa code |30( for calculations 
of elastic electron- nucleus scattering. The predictions ob- 
tained from the charge densities computed with various 
selected Skyrme forces and modern relativistic mean field 
parameter sets have been compared with existing exper- 
imental data about elastic electron scattering in several 
stable nuclei. 

A suitable quantitative comparison among the theo- 
retical predictions is established by introducing, for each 
nucleus and beam energy, a normalized square weighted 
difference. While all the considered effective interactions 
describe qualitatively well the experimental DCS data, 
the quantitative analysis in terms of the with respect 
to experiment shows some differences among the various 
theoretical DCS, especially for large scattering angles. 
These differences are related mainly with the different 
behavior of the mean field densities in the inner region 
of the nuclei. For the investigated nuclei and energies, 
one always finds a theoretical charge distribution whose 
is at least similar, and in some cases better, than the 
one obtained with the experimental charge density that 
has been fitted to reproduce the measured DCS data. 
Among the nuclear interactions analyzed, the DCS val- 
ues calculated with the charge densities obtained from 
the covariant DD-ME2 and FSUGold parameter sets are 
the ones that tend to give overall better agreement with 
the experimental values considered. 

A more challenging test on theory is provided by the 
analysis of the experimental data about relative differ- 
ences of DCS between pairs of neighbor nuclei. The cal- 
culations based on the mean field charge densities reason- 
ably follow the global trends shown up by the experimen- 
tal measurements. The quantitative analysis, however, as 
in the case of the previously investigated DCS, highlights 
some deficiencies of the mean field densities in their inner 
region. 

We have used the mean field charge densities obtained 
with the relativistic G2 parameterization, by the reasons 
discussed previously, as the input baseline to study the 
elastic electron-nucleus scattering along the tin and cal- 
cium isotopic chains. Calculations have been performed 
from the proton to the neutron drip lines. Our aim has 
been to extract general trends, according to current mean 
field theories of nuclear structure, about the behavior 
that may be expected from real electron-nucleus scatter- 
ing experiments in exotic nuclei, in the low-momentum 
transfer region. Such experiments are envisaged at FAIR 
[13 and the SCRIT [11, [11] project in the nearby fu- 
ture, using exotic nuclei provided by radioactive isotope 
beams. We have confined our study to medium and heavy 
mass isotopes where the mean field approach can be ap- 
plied. Surely, the experimental program will investigate 



not only heavy exotic nuclei but it will address also scat- 
tering from light exotic nuclei. The theoretical investi- 
gation of these light nuclei, however, demands a sophis- 
ticated microscopic treatment of scattering to deal with 
the underlying shell model structure and the possible oc- 
currence of halos, which is beyond the methodology of 
the mean field study of the present work. 

First, we have computed for each isotope of the inves- 
tigated chains the squared electric charge form factor. It 
has been obtained as the ratio between the DWBA DCS 
calculated with elsepa and the Mott DCS. We have 
checked that |i^DWBA(Q)P defined in this way is relatively 
independent of the energy of the electron beam up to mo- 
mentum transfers q « 1-1.5 fm"-"^, even for nuclei as large 
as ^"*Pb. Second, we have fitted the mean field charge 
densities by two-parameter Helm distributions. In doing 
so, we have adjusted the effective momentum transfer 
correction to reproduce the mean field rms charge radii, 
on average, along the isotopic chain. We have made an a 
posteriori check that a DWBA calculation of the elastic 
electron scattering using as input the equivalent Helm 
charge density, is in excellent agreement with the results 
computed with the original mean field charge density up 
to momentum transfers q w 1-1.5 fm~^. 

We have paid special attention to the value of the 
square of the DWBA electric charge form factor at the 
momentum transfer where its first minimum, in medium- 
mass nuclei, or its first inflection point, in heavier nu- 
clei, appears. We have studied how it evolves along 
the isotopic chains of tin and calcium. Interesting lin- 
ear correlations between the value of the momentum 
transfer at the first minimum (or IP) of |i^DWBA('?)P 
with the mass number of the isotopes of the chain, with 
the effective momentum transfer at the first minimum 
of |Fpwba(9)P, and with the rms radius of the charge 
distribution, have been discussed. Also, a linear corre- 
lation between |Fdwba(9)P and a'^q^ computed at the 
first minimum or IP, where a is the Helm parameter that 
accounts for the surface thickness of the nuclear density, 
has been found in the studied chains. 

The analysis described in the present paper could po- 
tentially be useful for future electron-nucleus elastic scat- 
tering experiments. If the experimental data are available 
for two or more isotopes of a given chain, the aforemen- 
tioned linear correlations would provide, for an unknown 
nucleus of the chain, a hint on the value expected for the 
square of the experimental electric charge form factor at 
its first minimum, and for the momentum transfer where 
the latter occurs. The parameters of the Helm charge 
density distribution of the unknown isotope could be es- 
timated by means of correlations such as those displayed 
in Figs. M and [H with the help of Eqs. ^ and (fT5|l . 
Also, if the value of the squared modulus of the form 
factor is determined experimentally at its first minimum, 
the charge density in the Helm model can be sketched 
from similar correlations to Figs. [51 and [T^ together with 
the correlation of the type depicted in Figs. [D and [T31 
The use of more elaborated versions of the Helm model 
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[stI . [6^ . [69| that take into account the central depres- 
sion of the charge density, should allow one to extend 
the domain of validity of our method up to larger values 
of the momentum transfer. Work in this direction will 
be undertaken. 
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